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ABSTRACT 

One of the key open questions in cosmology today pertains to understanding when, where and 
how super massive black holes form, while it is clear that mergers likely play a significant role in 
the growth cycles of black holes, how supermassive black holes form, and how galaxies grow around 
them. Here, we present Hubble Space Telescope WFC3/IR grism observations of a clumpy galaxy at 
z = 1.35, with evidence for 10^ — IO^Mq rapidly growing black holes in separate sub-components of 
the host galaxy. These black holes could have been brought into close proximity as a consequence of 
a rare multiple galaxy merger or they could have formed in situ. Such holes would eventually merge 
into a central black hole as the stellar clumps/components presumably coalesce to form a galaxy 
bulge. If we are witnessing the in-situ formation of multiple black holes, their properties can inform 
seed formation models and raise the possibility that massive black holes can continue to emerge in 
star-forming galaxies as late as z = 1.35 (4.8 Gyr after the Big Bang). 

Subject headings: galaxies: Seyfert; galaxies: high-redshift; galaxies: active 



1. INTRODUCTION 

Where supermassive black holes come from and under 
what conditions do they form are some of the currently 
open questions in astrophysics. Observational data that 
directly addresses this questi on at the earliest tim es are 
just becoming available (e.g.. iTreister et aI]|201lD . How- 
ever, it is unclear whether supermassive black holes con- 
tinue to be born throughout cosmic history. 

Most theoretical models for the formation of seed black 
holes in the Universe have focussed on early formation 
prior to the enrichment of the Universe by metals. In this 
Letter, we present the discovery of a triple active galactic 
nucleus (AGN) system in a galaxy at z = 1.35 featuring 
three fairly massive, efficiently-accreting black holes each 
hosted in a sub-component or clump of the host galaxy 
(Figure [T|). These three black holes could plausibly have 
formed either during early epochs, growing via merger- 
induced accretion, or more recently, growing in situ from 
seed black holes that collapsed within the preceding few 
hundred million years. The in-situ formation scenario 
is likely due to the expected rarity of the triple merger 
alternative, but it also has more wide-ranging implica- 
tions. We describe this unusual Triple AGN and the 
possible implication for black hole growth in this Letter. 
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Fig. 1. — GOODS Viz composite image {V - blue; I— green; z 
- red). The morphology is irregular and clumpy. There is no sign 
of any tidal tails indicative of a recent interaction or merger even 
at the imaging depth of the GOODS data. 



The discovery presented in this Letter was enabled 
by the new infrared grism capability of the Wide Field 
Camera 3 (WFC3) on board the Hubble Space Telescope 
(HST) which allows spatially resolved spectroscopy in the 
rest-frame optical of faint, moderate-redshift galaxies. 
Throughout this Letter, we assume a ACDM cosmology 
with hf) = 0.7, = 0.27 and Q\ = 0.7 3, in agreement 
with recent cosmological observations (jHinshaw et al.l 
[2001 . 



2.1. 



2. DATA & ANALYSIS 
Discovery and Basic Properties 



We discovered the galaxy in the Chandra Deep Field 
- North (CDFN) while examining t he WFC3/IR gr ism 
data of Chandra X-ray sources. iBarger et al.l (|2008[) re- 
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Fig. 2. — F140W image (left) and G141 grism dispersed image {right). In the F140W image, we label the four components of the galaxy 
A, B, C and D, respectively. The dispersed image shows emission lines. In the dispersion direction, components B and C are blended, but 
since B is brighter, the spectrum is dominated by B (see Figure |3)l. 
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Fig. 3. — Extracted 1-D spectra of components A, B (+C) and 
D. In each panel, we label the main emission lines. The large 
[Olll]/H/3 ratio, together with the [Sll] line identifies components 
A, B and D as AGN. The fits to the emission lines are shown in 
red. 



port a spectroscopic redshift oi z = 1.3550 (ID 1163), 
with ground-based coordinates of (J2000) 12:36:52.77 
+62:13:54.8. It is listed as an individuall y detected X-ray 
source (ID 272) in the 2 Ms catalogue of lAlexander et all 
(|2003), who report an X-ray flux of ^o.5-8 keV = 1-3 x 
2^0-16 j,j.g g-i g^jj^ hardness ratio of HR = 

{H-S)/{H + S) 0.03. At the redshift of the source, 

this flux corresponds to a modest observed luminosity 
of io.5_8 kcv = 1-4 X 10^^ erg s~^ and an obscuring 
colum n density of A^h ~ 5 x 10^^ cm~^ (jTreister et al.l 
|2009| ). which corresponds to a obscuration conversion of 
1.5 for the observed 0.5-8 keV band (rest-frame 1.2-19 
keV) for a typical AGN spectrum. The (optical) galaxy 
is larger than the Chandra point spread function but the 
low count number does not allow us to associate counts 
to individual components. The Hubble Space Telescope 
optical image is shown in Figure[l]and the infrared image 
and grism data are shown in Figure [2] 



We retrieved the WFC3 imaging and grism data of 
the CDFN field^ from the archive and processed them 
with standard software. The detection F140W image 
was reduced using the ST SDAS pyraf task multidrizle 
(|Koekemoer et all I2002D . The dispersed (G141 grism) 
data wer e reduced using the aXe slitless spectroscopy 
package (jKiimmel et al.lf2009( ). 

In Figure [2] we show both the undispersed F140W {H- 
band) image as well as the dispersed G141 grism image. 
The source shows two complexes of two components each, 
which we label A, B, C and D. The orientation angle of 
the dispersed observation means that while two of these 
components, A and D, are cleanly dispersed, two others, 
B and C, contaminate each other. Fortunately, source B 
is significantly brighter than C, so the dispersed data are 
dominated by component B. We extract three different 
spectra from the dispersed image containing source A, 
B(-f C) and D, respectively, and analyze them separately 
(Figure El). 

For each of the three extracted spectra, we measure 
the [Oiii] A5007 and II/3 line strengths, as the [Oiii]/II^ 
ratio is a diagnostic of the domin ant source of ioniza- 
tion (AGN vs. s tar formation; [Baldwin et al.l 119811 : 
IKewlev et al1l200l . At the resolution of the WFC3 G141 
grism, the Ha line is completely blended with the [Nil] 
A6548 and A6583 doublet surrounding Ha. For AGN 
photoionization, the Ha and [N Ii] doublet become com- 
parable in strength, so we cannot use the blended line as 
a diagnostic as the [Nil] /Ha ratio. We account for the 
blending of [O ill] A4959 and A5007. 

We therefore measure only the [O iii]/H/3 ratio for each 
of the three components by fitting Gaussians using the 
mpf it algorithm (Markwardt 2009) from the extracted 
grism spectra in Figure[3] All three spectra show a strong 
[Qui] A5007 line and a weak H/3 line yielding line ra- 
tios of logio([0 iii]/H/3)~ 0.4 — 0.5 which are commonly 
associated with photoionization by AGN. In two of the 
three components (A and B), the [Sii] AA6717, 6731 line 
is clearly visible supporting the identification of these 
components as AGN and ruling out the low-metallicity 
starburst scenario. The spectrum of Component D is too 
noisy to reliably detect the [Sii] line but there appears 
to be a marginal detection. We report these results in 
Table [H 

Theoretical models could produce such a high 
[Qui] /H/3 ratio by an extreme combination of signifi- 
cantly sub-sol ar metallicity and v ery high ionization pa- 
rameters fe.g.. IKewlev et al.ll200lD but this alternative is 
not consistent with the X-ray data. Specifically, if we 
attribute aU the flux in the blended Ha+[Nii] A6584 
region to a single Ha line to star-formation, the im- 



2.2. WFC3/IR Imaging and Grism Data 



* Taken as part of the Cycle 17 Proposal ID #11600. 
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TABLE 1 

Component Properties 



Component 


L[Oni] ASOO?" 
erg s^-"- 


H/3 
erg s ■'- 


logioCIOiiil/H/?)" 


Stellar Mass 
Mq 


Black Hole Mass 

Me 


A 
B 
C 
D 


2.7(±0.14)xl0*2 
2.3(±0.14)xl0*^ 

6.8{±1.3) xlO"'^ 


6.8(±1.2)xl0''i 
6.6(±1.0)xl0-*i 

1.75(±0.8)xl0''^ 


0.53 ± 0.1 
0.44 ± 0.1 

0.37 ± 0.2 


9.6t"77 X IQS 

2.8lo48 X 10^ 
SAtfl X IQS 


3.1x10'' 
1.2x10^^ 
(2.9x10'*)'' 
1.0x10^ 



The [Olll] A5007 line is corrected for the blended [Olll] A4959 component. 
^ The grism data do not allow determining whether component C contains an AGN. 



plied maximum star formation rates are 11.0, 6.9 and 
2.1 Mp, yr~i for components A, B and D, respectively 
ennicutli 119981 ). while the star formation rate implied 
by th e X-ray luminosity is 420 Mq yr"^ (|Ranalli et alj 
|2003[ ). Furthermore, the star- formation rates and stel- 
lar masses of each clump imply gas-phase metallicities 
that are too high to explain the observed [Oiii]/H/3 ra- 
tios were they due t o star formation. Follow ing the 'fun- 
damental plane' of IMannucci et al.l (|2010[ ) and assum- 
ing that all the blended Ha-|-[Nii] A6584 is due to star 
formation, we compute 12-l-O/H for the three compo- 
nents of 8.49, 8.77 an d 8.82, respectively. According to 
IMaiolino et al.l (j2008D . these high gas phase metallicities 
would imply logio[0 iii]/H^ ratios of 0.41, 0.03 and -0.06, 
i. e. high metallicity. The first of these is marginally be- 
low the measured logio[0 iii]/H/3 — 0.53 ± 0.1 for com- 
ponent A, but the other two are completely inconsistent. 
As with the other arguments here, a low metallicity star- 
burst cannot be completely ruled out but the preponder- 
ance of the evidence points to three AGN. Future obser- 
vations, e.g. with JWST, could settle the debate. 

2.3. ACS Imaging and Morphological Decomposition 

We retrieved the ACS imaging data in the B, V, I and 
z bands taken as part of the GOODS survey from the 
archive and determined the morphologies and luminosi- 
ties of each component source using th e two-dimensiona l 
image decomposition program GALFIT (jPeng et al.ll2002[ ). 
We find that each of the components in all ACS bands 
have morphologies consistent with disks. Based on anal- 
ysis of more than 50, O OP simulated AGN-|-ho st galaxy 
morphologies at z < 1 (Simmons, fc Urrvl[2008l ). this in- 
dicates that no more than 20% of each object's rest-frame 
UV light comes from a bulge. 

In order to investigate whether any of the components 
have detected central point sources in the ACS images, 
we fit each component with both a single Sersic com- 
ponent and a Sersic -I- point source. We assess the ro- 
bustness of the point-source detection via analysis of the 
goodness-of-fit parameter and examination of the fit 
residuals. Only for component D is the improvement 
in marginally significant, and the fit is unable to con- 
verge to physically realistic parameters without a central 
point source. We therefore determine that component D 
has a detected point source in the rest-frame UV bands, 
supporting the cl aim that D is an AGN. F rom our exten- 
sive simulations ([Simmons fc Urrvl [20081) . the likelihood 
of a falsely detected point source is very small (< 1%), 
while there is a 5% chance of missing a central point 
source. While the addition of a central point source does 



not improve the fit of component A, the addition of a 
second Sersic parameter does. Component A is best fit 
in all bands by two disks, one compact and one more 
extended. 

2.4. Black Hole Masses 

Locally, black hole mass is correlated with galaxy 
bulge velocity dispe rsion ([Ferrarese fc Merr"itt[ [20001: 
Gebhardt et al.l |20C)0D and stellar mass ([Haring &: Ri^ 
20041 ). We obtain a photometric spectral energy distri- 
bution (SED) of all four clumps using the ACS BVIz 
and WFC3 F140W images and fit a two-component star- 
forma tion history to each cl ump to estimate its stellar 
mass ([Schawinski et alll2007D . From these stellar masses, 
reported in Table [H we infer b l ack h ole masses using the 
local relation of H iring fc Rixl ([20041 ) , yielding masses in 
the range of 3 x 10'^ to 2 X lO'^M©. We caution that there 
are two major systematics that could affect these mass 
estimates: (1) the mass in each clump is assumed to be 
all bulge, and (2) the local black hole - bulge r elation 
may evolve to z = 1.3 (e.g. iBennert et al.l[2011[ ). The 
first consideration means that these masses are upper 
limits to the true black hole mass; since the morpholog- 
ical analysis indicates that any bulge is less than 20% of 
the observed light, the black hole masses are plausibly at 
least 5 times smaller. 

2.5. Eddington Ratios 

The [O III] A5007 luminosity ha s been used as a prox y 
for black hole accretion rate ('e.g.. iHeckman et al.ll200'i) . 
We correct the observed [Oiii] luminosity measured in 
the grism spectra for the contribution of the blended 
[Oiii] A4959 contribution to use it as a basis for calcu- 
lating the accretion rate and therefore Eddington ratio 
of each black hole. We scale L[Oiii] to an X-ray lumi- 
nosity using the empirical correlation between the two 
for moderate redshift AGN (a factor of ~ 10; Figure 5 
of Tr ouille fc Ba rger 2 01Qi) . We then scale this X-ray 
luminos ity to a bolometric l uminosity. The local rela- 
tion by iMarconi et al.l ([2004f ) suggests a factor of ^ 20 
which is supported by r ecent measurements of z ~ 1 
AGN by[S immons et al.l ([20111 ). Using bolometric lumi- 
nosities obtained this way results in Eddington ratios of 
L/L^dd — 1.4, 0.3 and 0.09 for components A, B and 
D, respectively. We note that these Eddington ratios 
are likely underestimates due to the assumption that the 
whole stellar mass in each component represents pure 
bulge light and the 20% limit on the bulge light would 
lower the black hole masses, and raise the Eddington ra- 
tions, by a factor of 5 respectively. 
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Fig. 4. — Growth rate of black hole seeds as a function of time 
assuming Eddington-limited accretion, a radiative efficiency of e = 
0.1 and tEdd = 0-45 Gyr. The shaded gray area indicates the range 
of black hole masses estimated here for the Triple AGN. This mass 
range can be reached comfortably via Eddington-limited growth 
within ~few hundred million years for a large range of seed masses, 
and under 100 Myr for massive seeds. 

The observed Chandra X-ray properties are derived 
from 11 photons which cannot be locahzed to individ- 
ual clumps and so represent the summed emission from 
the three AGN. It is possible that the unobscured X-ray 
light is dominated by the AGN in component D which 
features a point source, with the other two AGN being 
more obscured. This yields an Eddington ratio of 0.03 
for component D, assuming a bolometric correction fac- 
tor of 20. The AGN in components A and B are likely 
more obscured but still visible in [O ill] and therefore 
intrinsically more luminous. 

3. DISCUSSION 

We have described a galaxy at z = 1.35 consisting 
of four individual sub-components visible in both ACS 
BVIz and WFC3/IR F14GW images. From slitless spec- 
troscopy with WFC3/IR in grism mode, we identified ac- 
creting black holes in three of these sub-components. We 
estimated black hole masses for each clump and found 
low black hole masses in the range of 10^ — 10'' or 
even smaller. 

3.1. Origin of Seed Black Holes and Rapid In-Situ 
Growth? 

The in-situ formation scenario is plausible due to the 
fact that these black holes are low mass and are accret- 
ing at a high Eddington rate. As we estimate below, 
even a low mass black hole seed can reach the observed 
mass range in a short time without requiring any extreme 
super-Eddington accretion episodes. To illustrate, we es- 
timate below the growth as a function of time for a black 
hole accreting at the Eddington limit : 

M(t) = M(0)expfi^-^V (1) 

where the radiative efficiency e w 0.1 and the characteris- 
tic time scale for Eddington-limited growth is t-^dd = 0.45 
Gyr. The growth curves for black hole seeds ranging 



from 1 — IO^A/q in Figure H] show that even very low 
mass seeds of ^ 10^ Mq can comfortably reach the ob- 
served black hole mass range of 10^ — lO^M0in ~500 
Myr, while a more massive ^ IO^Mq initial seed will 
reach the observed range in under 100 Myr. Their lo- 
cation at the centers of clumps embedded in a larger 
galaxy means that there likely is gas to sustain accretion 
to the observed masses. Similarly, the growth times for 
massive seed models approach the dynamical timescale 
of the host galaxy. 

Various mechanisms for producing black hole seeds, es- 
pecially massive black hole seeds at extremely high red- 
shift to account for z ^ 6 quasars , have bee n proposed 
(iLodato fc NataraianI 120061 120071: see also iNataraianI 
l2011f ). Two general channels have been theoretically 
proposed and extensively studied: (1) the remnant 
of a Population III star, (2) massive primordial gas 
disks that can either collapse directly into a black hole 
or fragment into a dense star cluster that collapses 
into a black hole via runaway collisions and s tellar- 
dynamical processes (iLo dato & Nataraia;^ 120061 120071 : 
iDevecchi fc Volonterill2009t iVolonteril 12010 1. The loca- 
tion of the three AGN at the centers of clumps embedded 
within a galaxy does not allow us to distinguish between 
these scenarios. 

If indeed these observed black holes formed in situ as 
recently a.s z = 1.35, our theoretical models should be 
re-examined, as all the channels considered so far have 
been worked out in the context of primordial gas disks. 
The presence of metals at late times will rapidly enhance 
cooling and the outcomes for this case are yet to be ex- 
plored. Even with th e one proposal to fo r m bla ck holes 
as late as z = 3 — 4, iJimenez &: Haim"anl ()2006[ ) require 
pockets of unenriched gas to exist, but it is unlikely that 
such conditions would survive to z = 1.35. A channel 
allowing seed formation at non-primordial metallicities 
may be required. 

We note that the possibility that these black holes 
formed significantly earlier, perhaps even via the Pop III 
mechanism, cannot be ruled out. In this case, the present 
accretion event represents a late-stage re-activation. 

3.2. Three- Way Merger or Clumpy Galaxy? 

With the current data, we cannot fully determine 
whether the host galaxy of the Triple AGN is a ma- 
jor merger of three individual progenitor galaxies, or 
whether it is a clumpy galaxy. Clumpy galaxies 
have been observed iii deep Hubble observatioii s (e.g . 
Moustakas et al.| 120041: lElmegreen fc Elmegreeiil 120051: 
Elmegreen et al.ll2009f) while mergers of multiple galaxies 
are rare, at least in the local Universe (two orders of mag- 
nitude less frequent than binary mergers, see Darg et aL| 
MvH ). Neither the ACS optical nor the WFC3 IR im- 
age show clear signs of tidal tails, merger debris, or close 
companions so even with the higher frequency of mergers 
at high redshift, so it seems more likely that the Triple 
AGN is a clumpy galaxy. 

There are only two cas es of Triple AGN sys tems 
reported in the literature (jDjorgovski et al.l 120071 and 
iLiu et all [20^1 11. These multiple AGN systems are iden- 
tified as merging systems. Each progenitor is thought 
to have brought in an already existing black hole and 
the double AGN represents the stage of the merger prior 
to the final coalescence into a single central black hole. 
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The properties of both these cases are different from the 
object presented here: the z ~ Triple AGN contains 
massive > IO^Mq black holes radiating at significantly 
sub-Eddington rates, while the z ~ triple quasar has 
high black hole masses and high luminosity. 

Clumpy or irregular galaxies are seen in deep Hub- 
ble Space Telescope fields, rise in abundance with red- 
shift, and make up > 50% of ga laxies at z > 1 
(iMoustakas et alj 12004 lElmegreen fc Elmcgree^ 120051: 
IConselice et alJ l2005t lElmegreen et alJ 120091) while the 
fraction of merg i ng galaxies rema ins below ~ 10% 
(jBell et alJ [200a iLotz et all [2008I ). The clumps in 
this common type of star-forming galaxy are thought 
to be fragmenting gas disks where individual re- 
gions become self-gravitating, star-forming clumps (e.g. 
NoEuchi 1999; Immcli et al. 2004; Bournaud et al. 2007; 
Elmegrecn et al.ll2008al |bl: IBournaud et al.ll2011l) . 

Spectroscopic studies of high redshift star-forming 
galaxies reveal that many of them are turbulent, gas-rich 
disks with clumps rapidly forming without the need 
for mergers (iGe nzel et al. 2008; Fo rster Schreiber et al.l 
[20091 IDaddi et all 120101: iForster Schreiber et al.l 120111 ). 
Theorists have reproduced such gas-rich disk galax- 
ies with the clumps emerging by fragment ation as 
sub-clu mps become self-g ra vitating (e.g. Noguchil 
19991: llmmeh et al.l |2004t iBournaud et al.l 120071: 
Elmegreen et al.l l2008allbl : IBournaud et a l.' '20li| ). If 
the clum ps are seeded with black holes. [Elmegrec n et al.l 
(12008b) have shown that as the stellar clumps merge 
into a bulge, the black holes similarly coalesced into 
a supermassive black hole at the center, preserving 
the local M - a re lation (jFerrarese fc MerrittI 120001: 
IGebhardt et"aI1[2000l) . 

As dense, gas-rich, collapsing regions of space, they are 
viable birth places for the direct formation of supermas- 
sive black holes and our discovery of the z — 1.35 Triple 
AGN shows that - if the host really is a clumpy galaxy - 
they can host growing black holes inside their constituent 
clumps. 

4. SUMMARY 

We have discovered a triple source a.t z — 1.35 which 
evidence suggests are three black holes, having masses 



10^ — 10*" Mq or smaller, growing at a significant frac- 
tion of their Eddington luminosities. These three black 
holes were possibly born in their current locations in sub- 
galactic clumps. This raises the possibility that super- 
massive black holes were not only formed in the very 
early Universe, but continue to form within galaxies as 
late as z = 1.35 (4.8 Gyr after the Big Bang) where 
the extreme low metallicity invoked in some early Uni- 
verse black hole seed formation models are unlikely to be 
present. While the present observations of the system do 
not confirm or reject any specific seed-formation scenar- 
ios, its discovery raises for the first time the possibility 
of observing seed formation in action. 

The Triple AGN presented here was discovered 
serendipitously and is barely detected in the X-rays in 2 
Ms Chandra observations; only slightly less massive black 
holes may thus lurk just below the detection limit of deep 
X-ray AGN surveys. While this object may represent a 
pathological case, more systematic searches for low-mass 
black holes in z ~ 1.3 galaxies will reveal whether this 
is a common mode of black hole formation and growth 
in the Universe and the physical conditions within the 
clumps may shed light on models for black hole seed for- 
mation. If these are indeed low mass black holes that 
formed in situ, we require new theoretical models for the 
late formation of black holes in the Universe. 
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